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Abstract 

We calculate the Mellin moments of the 0{a1) coefficient functions for the unpolarized 
and polarized fragmentation functions. They can be expressed in terms of multiple finite 
harmonic sums of maximal weight w = 4. Using algebraic and structural relations between 
harmonic sums one finds that besides the single harmonic sums only three basic sums and 
their derivatives w.r.t. the summation index contribute. The Mellin moments are analyti- 
cally continued to complex values of the Mellin variable. This representation significantly 
reduces the large complexity being present in x-space calculations and allow very compact 
and fast numerical implementations. 



1 Introduction 

The production of hadrons in annihilation processes such as e^e~ or quark-antiquark and gluon- 
gluon scattering is described by parton-hadron fragmentation in the final state of these pro- 
cesses. Due to the factorization theorems perturbative Quantum Chromodynamics describes 
the hard scattering sub-processes at the level of leading twist, while hadron formation requires 
non-perturbative functions, the fragmentation densities Dp{z,M'^). At leading order in the 
coupling constant they describe the probability that a final state parton p involved in the hard 
scattering process produces a hadron H. Here z denotes the momentum fraction of the parton 
in the hadron and M^ is the respective factorization scale. At higher orders the fragmentation 
process is described by fragmentation functions which are formed of Mellin convolutions of the 
fragmentation densities and the corresponding time-like Wilson coefficients. 

The process of e"^e~ annihilation into hadrons provides the cleanest environment to measure 
the fragmentation functions and to extract the individual fragmentation densities : 

e+ + e- ^H + X . (1) 

By crossing from s- to t-channel this process transforms into deeply inelastic scattering of an 
electron off the hadron H : 

e- + H^e-+X. (2) 

Both processes are related by this transformation and the question arises for similarities or 
transformation properties of characteristic quantities contributing to each of the reactions. This 
question has been raised early [1] concerning both the non-perturbative and perturbative compo- 
nents. Crossing of the non-perturbative fragmentation- and structure functions is not expected 
to hold [2]. However, one may ask the question for the perturbative part order by order in 
the strong coupling constant a^ = as/{4:7i). Here it is not expected that the crossing holds 
for the Wilson coefficients and splitting functions defined in a specific factorization scheme, as 
e.g. the MS scheme, individually. However, physical quantities as physical evolution kernels Kij 
obey the crossing relations. These functions describe the scale evolution of observables and are 
scheme-independent. To NLO they were derived in [3-6] for polarized and unpolarized targets 
(final state hadrons) referring to different observables as the structure (fragmentation) functions 
and their slope and combining the longitudinal and transverse projections F2 and F^. Recently 
the scheme-invariant evolution kernels were derived at NNLO in [7] for the space-like case^. In 
Ref. [6] the conditions were derived up to the 2-loop Wilson coefficients. The crossing of the 
2-loop splitting functions was studied [8] based on Refs. [9]. 

To be capable to construct physical evolution kernels in the space- and time-like cases, the 
complexity of the respective integrals requires to choose the Mellin space representation. For 
the whole set of the time-like Wilson coefficients so far only the x-space representation has 
been available [10-13]^. In this paper we derive the Mellin-representation of the 2-loop Wilson 
coefficients. As observed in case of other single-scale quantities at the 2- and 3-loop level before 
[14-18] a considerable structural simplification for these quantities is obtained w.r.t. the number 
of contributing functions. In the Mellin-space representation the evolution equations can be 
solved analytically and fast and precise numerical implementations are obtained [19,20]. The 2- 
loop Wilson coefficients depend only on a few universal basic functions, for which representations 
in the complex Mellin-plane have to be derived. 



^The structure in the time-like case is analogous. 

^In the present analysis typos contained in the print were accounted for. 



The paper is organized as follows. In section 2 the differential cross sections for the process 
e^e^ ^ H + X are summarized and expressed in terms of fragmentation densities and coefficient 
functions to 0{a1) in the unpolarized and the polarized case. In section 3 we derive the Mellin 
moments for the coefficient functions and discuss their structure w.r.t. the set of basic functions 
emerging. It turns out that this is a subset of the functions needed in the space-like and related 
cases [16, 18]. The analytic continuations for these functions to complex values of the Mellin 
variable A^ was given before in [21]. Section 4 contains the conclusions. In the Appendices A 
and B the explicit representations of the Wilson coefficients in the unpolarized and polarized 
case in Mellin space are presented. 

2 Coefficient Functions 

Let us consider the differential scattering cross section for e^e~-annihilation into a hadron H 
and the hadronic remainder part X 

e+ + e" ^ 7, Z -^ if + 'X' . (3) 

This reaction is either studied for unpolarized or polarized leptons. If in the latter case the 
polarization of the produced hadron H is measured, one can form the polarization asymmetry. 
In the following we will consider both the case of unpolarized scattering and the final state 
polarization asymmetry in the polarized case. 

2.1 Unpolarized Case 

—- = - (l + cos^^) — ^ + -sm2^— ^ + -cos^— -^ . (4) 

dxd cos 8 ax 4 ax 4 ax 

Here x = 2p.q/Q'^ denotes the scaling variable with Q^ = q^ > and 6 is scattering angel of the 
produced hadron. p denotes the momentum of the produced hadron and q is the virtuality of the 
exchanged 'y{Z). x denotes the fraction of beam energy carried away by the produced hadron H. 
The differential scattering cross section (4) consists of three contributions o"^^ ^ corresponding 
to the transverse (T), longitudinal (L) and the asymmetric (A) part due to 7 — Z interference 
and pure Z-boson exchange. The individual contributions to the differential cross sections are 
given by 
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=1 
with k = T,L. The asymmetric cross section is given by 

^ = j:Ar\Q^)Dl{x,M^)^CzU^) . (6) 

1=1 ^ ^ 

Nf denotes the number of flavors. The Mellin convolution ® is deflned by 

[Ai^B](x)= / dxi / dx25{x — xiX2)A(xi)B(x2) ■ (7) 



In (5, 6) the flavor singlet and non-singlet combinations of the non-perturbative fragmentation 
densities Di[z,M'^) contribute, which are given by 



Di{z,M') 



f 1=1 



D^sAz,M' 



= Df{z,M')+Dj{z,M')-Ds{z,M^), 
DUz^M') = Dl'iz,M')-Djiz,M'), 



(9) 
(10) 



with M^ the factorization scale. The electro-weak pointlike cross sections are 
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(12) 



Z{Q^) =Q^-Ml + iMzTz 



(13) 



and Mz, T z denote the mass and width of the Z-boson. The electric charges are cl = —1, e„ = 
2/3, Cd = —1/3 and Cv,{A),l, Cv,{A},i are the electro-weak couplings of the charged lepton and 
the quarks, respectively. Nc = 3 denotes the number of colors. Correspondingly the asymmetry 
factor reads 
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The non-singlet and singlet coefficient functions are denoted by C? / Jx, Q^/M^) and 



^ll{K)i.^^Q'^lM'^)- To 0{al) they are given by 
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Here a^ = Q;s/(47r) = g'l/{ATiY and (^^ denotes the strong coupling constant. P-^ are the the 
/ + 1-st order QCD sphtting functions and /^o = 11Ca/3 - (4/3)TRiV/ with Ca = S,Tr = 1/2 for 
SU{3)c- The superscripts 'ici' and 'nid' refer to identical quark contributions resp. the remainder 
part. The corresponding Wilson coefficients to 0{a'^) were calculated in Refs. [10-12]. 

2.2 Polarized Case 



The polarization asymmetry for the fragmentation process was calculated in [13] to 0(0^). The 
hadronic tensor of the process is obtained from the hadronic tensor in polarized electro-weak 
deeply inelastic scattering [22] by crossing from t to s-channel. In general five fragmentation 
functions gf{x, Q^)\^=i contribute to the scattering cross section. If we limit the analysis to the 
level of twist-2 fragmentation functions the cross section is determined by (7^45. Here g^^^ are 
non-singlet fragmentation functions which contribute to the 7 — Z and | Z | ^-exchange terms. 
The corresponding coefficient functions are the same as for the fragmentation functions F2 and 
Fi, respectively, in the unpolarized case due to the structure of the hadronic tensor, cf. [22]. In 
the following, we will therefore consider the case of pure photon exchange only. If the scattering 
cross section is integrated over the azimuthal angel of the produced hadron the differential 
cross section is given by 



dxd cos 6 dxd cos 6 
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The fragmentation function g^ is given in terms of polarized fragmentation densities A 
and Wilson coefficients AC;^ ^ , 
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The polarized fragmentation densities are 
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The respective singlet and non-singlet combinations are given by 
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The polarized non-singlet and singlet coefficient functions ACf ^(x, Q"^ /M"^) and 
AC^^g(^)(a;, QV^^) to 0{a'^) are given by 
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Here AP/ ■ denote the polarized splitting functions. The non-singlet Wilson coefficients for the 



«j 



structure functions (71,4,5 are the same as in the unpolarized case for the corresponding tensor 
structures, i.e. 

(32) 
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The structure functions (745 are pure non-singlet. 



3 Mellin Moments 

The coefficient functions to Ola"^^) in x-space can be expressed in terms of Nielsen integrals 

SnAy)= \ \„ , / -ln-\z)ln^{l -zy), (35) 

(n- ly.pl Jo z 

partly with y being a rational function of x [10-13]. Here a relatively large number of functions 
contributes. The weight w of these functions is defined by w = p + n, where any power of 
a logarithm counts for w = 1 as is the case for the denominators 1/x, 1/(1 + x) and 1/(1 — 
x). In a product of functions the weights of the factors add. The Nielsen integrals cover the 
polylogarithms Sn~-\,\{y) = Li„(i/) and through dLi2{y) / d\ia{y) = — ln(l — y) the logarithm. 

For massless calculations it is suitable to consider the coefficient functions in Mellin space. 
The Mellin-transform of a function F[x) is defined by 

M[F](iV) = f dxx^-^F{x) . (36) 

Jo 

In particular, for the Mellin convolution [Fi (g) Fq\ {x) of two functions -Fi(x), F2{x) one obtains 

M[Fi ® F2]{N) = M[F,]{N) ■ M[F2]{N) . (37) 

In Mellin space fragmentation densities and the corresponding Wilson coefficients are connected 
multiplicatively. Mellin space representations have also the advantage that the QCD evolution 
equations can be solved analytically [19,20] which are numerically very precise and fast. 

The Mellin transforms of single-scale quantities in QCD can be represented in terms of nested 
multiple harmonic sums [23,30] and rational functions of the Mellin variable A^. In this way a 
standardization is obtained. Harmonic sums are described recursively by 

S,,a,,...,a,{N) = J2 ^^7^^a,...,a,(/) , (38) 

1=1 

where b,ai, . . . ,ak are positive or negative integers and A^ is a positive integer. Sa{N) is called 
single harmonic sum. 

The harmonic sums S'ai,...,afc(A^) form a Hopf algebra [24,25]. For each index set {ai, . . . , a^} 
one can determine the basis elements through which all harmonic sums belonging to this set can 
be expressed algebraically. The number of these basic sums is given by the number of Lyndon 
words of this set and can be calculated by Witt-formulae [25,26]. 

For the application of the Mellin-space representation to experimental data one needs to 
derive analytic continuations of the harmonic sums from integer values oi N, N - even or odd -, 
to complex values. While in the case of single harmonic sums these representations are known, 
they are more difficult to derive for multiple harmonic sums [17,21]. As finally the argument 
A^ will be complex one may consider further relations between harmonic sums for rational and 
real values of the argument A^ [27] beyond the algebraic relations. In this way the amount of 
basic sums is further reduced. In particular one may consider equivalence classes of functions, 
which are related by differentiation. If an analytic representation for a harmonic sum is found 
in the complex plane outside its singularities, any derivation is easily obtained. It is convenient 
to choose for the basic functions Mellin transforms of a suitable harmonic polylogarithm [28] or 
Nielsen integrals with argument x weighted by 1/(1 =F x). 



In former analyses of the two-loop anomalous dimensions, the two-loop Wilson coefficients 
for unpolarized and polarized deeply inelastic structure functions [18], the unpolarized and po- 
larized Drell-Yan process and scalar and pseudoscalar Higgs boson production [16], 0{a^) heavy 
flavor contributions to Fl{x^Q'^) [29] and the unpolarized 0(q;^) anomalous dimensions [17] the 
following minimal sets basic functions w.r.t the equivalence classes deflned above were found. 
We list the newly emerging basic functions. 

i) 1 loop anomalous dimensions and coefficient functions -.^ 
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ii) 2 loop anomalous dimensions : 
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iii) 2 loop coefficient functions : 
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In the case of the 3-loop anomalous dimensions further eight basic functions contribute [17]. 
The notation above follows Ref. [21]. Since 
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one obtains 



M[/(x)](iV)= m[\n\x)f{x)\{N), 
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Ar,(N) and A22{N) are represented by A^i^N) and Ais{N) in the sense of equivalence classes as 
outlined above. 

The Mellin transforms of the unpolarized and polarized time-like coefficient functions to 
O(a^) are given in Appendix A and B. Here we discuss their principal structure. Unlike the case 
for the 0(q;^) coefficient functions for deeply inelastic scattering, the process related through 
crossing from the time-like case considered in this paper to the space like-region [6], besides the 
single harmonic sums only three other basic function contribute. For the deep-inelastic Wilson 
coefficients [18] and the polarized and unpolarized Drell-Yan process and scalar and pseudoscalar 
Higgs-boson production cross sections in the heavy mass limit [16] five additional functions 
occurred. As the corrections to the latter hadronic processes are initial-state corrections they 
belong to the same class as space-like processes. Furthermore it is interesting to see, whether 

^In case of functions which diverge as a; ^ 1 the (...)+ prescription is understood. 



any of the harmonic sums contributing contains an index {—1} since those sums were generally 
absent in all cases analyzed before. Alternating single harmonic sums are given by 
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and iIj{x) = dln{r{x))/dx. In Appendix A and B the /3— function indeed occurs. However, it 
emerges always together with the function A^lN) : 
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in such a way, that it cancels. Therefore, all contributing harmonic sums do not contain the 
index { — 1} also in the present cases. 

The analytic continuations of the Wilson coefficients from integer values to complex values of 
N is now easily performed. The alternating single sums are expressed in terms of Eqs. (45,46). 
The non-alternating single sums obey the representation 



Sk{N) = (_l)(i+^)^^l-^^('=-i)(iV+l) + C(A;). 



(49) 
(50) 



The analytic continuations for the remaining functions A3{N),Ais{N) and A2i{N) were per- 
formed in [21] and A^^N) and A22{N) are obtained as derivatives (43,44). 

In the case of the longitudinal fragmentation functions a transformation to Mellin space was 
given in [31] recently using the transformation tables of one of the authors [30], see also [32]. In 
Appendix A and B we present the Mellin transform in such a way, that factors of (—1)^ still 
present in [31], Eq. (A5), do not occur, and use further compactifications compared to earlier 
results. As well known, the analytic continuation of the Mellin moments is carried out either from 
the even or odd integer values of N determined by the crossing relations q — > —q,p — *> p of the 
respective quantity, cf. [22]. Thus, factors of (—1)^ are not present in the analytic continuation. 

A FORTRAN-programme of the 2-loop Wilson coefficients for the representation of the frag- 
mentation functions for complex values of the Mellin variable suitable to extend analyses of the 
experimental fragmentation data [31,33] can be obtained by the authors. 

4 Conclusion 

The mathematical structure of the 2-loop Wilson coefficients for the polarized and unpolarized 
time-like fragmentation functions simplifies considerably in Mellin space expressing the result 
in terms of nested harmonic sums. Using algebraic and structural identities for the harmonic 
sums and their continuations to rational, real and finally complex values of the Mellin variable. 
The physical evolution kernels for deeply inelastic scattering and parton fragmentation into 
hadrons are related by crossing from t— to s— channel order by order in perturbation theory. 



9 



Since these quantities are rational functions of the respective Wilson coefficients and anomalous 
dimensions in Mellin space one expects a familiarity of the mathematical structure of the space- 
and time-like Wilson coefficients. While in the space-like case five non-trivial basic functions 
emerge besides the single harmonic sums, in the time-like case only three of them contribute 
to Ola"^). In x-space many more functions are required to express the Wilson coefficients. As 
observed for various other processes, harmonic sums with index {—1} do not contribute. The 
Mellin-space expressions allow to use the time-like Wilson coefficients in fast evolution codes for 
the experimental analysis of fragmentation function data. 

Note added. The time-like Wilson coefficients in the unpolarized case have been newly calculated 
in [34] very recently. 
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A Unpolarized Time— like Coefficient Functions 

In this appendix, we present the Melhn moments of unpolarized time-hke coefficient functions 
to Ofa?) : 



withF = NS,PS,S. 
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B Polarized Time— like Coefficient Functions 

In this appendix, we present the Melhn moments of the polarized time-hke coefficient functions 
to Ofa?). 



The individual contributions are 
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